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We formulate and approximately solve a specific many-body generalization of the Landau-Zener 
problem. Unlike with the single particle Landau-Zener problem, our system does not abide in the 
adiabatic ground state, even at very slow driving rates. The structure of the theory suggests that 
this finding reflects a more general phenomenon in the physics of adiabatically driven many particle 
systems. Our solution can be used to understand, for example, the behavior of two-level systems 
coupled to an electromagnetic field, as realized in cavity QED experiments. 
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The Landau-Zener (LZ) problem describes a paradig- 
matic situation in physics where two quantum levels cross 
each other in time. In its most basic form, the problem 
is represented by the Hamiltonian 

where t is time, g the coupling constant, and A the rate 
of change of the energy levels (here, as in the rest of the 
paper, we set h — 1). This Hamiltonian has two instan- 
taneous energy levels E± — ±^(\t) 2 + g 2 . Suppose in 
the distant past, t — > — oo, the system is in level E-. The 
goal then is to calculate the probability, P, to stay in E_ 
at t — > +oo. Solving the corresponding time dependent 
Schrodinger equation, Landau [T] and Zener [2_ found 

P = l-e-^, (2) 

as an exact answer to this question: if only the sweep- 
ing rate is slow enough, it is exponentially likely that 
the system will abide in its adiabatic ground state. For 
75 years, the Hamiltonian 0, and its solution ^ have 
been used to describe a huge spectrum of physical phe- 
nomena [3J. Subsequent generalizations of include an 
extension to a multi-channel environment wherein the 2 
by 2 matrix is replaced by a larger time-dependent ma- 
trix. However, common to all those problems [4] [5] [6] [7] 
is that only a finite number of degrees of freedom par- 
ticipate in the transition process (which manifests itself 
in transition probabilities of the same algebraic structure 
as in Eq. (g.) 

At the same time, there appears to be some interest in 
genuine [S] many-body generalizations of the LZ setup: 
fundamentally, one would like to know whether a slowly 
driven many body system will remain in its adiabatic 
ground state, in a manner resembling the single particle 
case Q . But there is also applied relevance to the gener- 
alization. A number of existing experimental setups pro- 
vide a perspective to actually probe the transition rates 



of a many-body LZ problem. Examples include systems 
of N two-level systems ('atoms', either real or artificial), 
coupled to a photon mode in a cavity [SJ HO] . In this 
case, time dependence might be introduced by changes 
of either the photon frequency (by changing the cavity's 
size), or the energy splitting of the two level systems (by 
applying a 'Zeeman' field.) Similar physics also arises in 
the context of polaritons, excitons coupled to a cavity 
electromagnetic mode Another phenomenon rele- 

vant to the present work is the observation of molecule 
production in an atomic gas experiment, due to sweeping 
through a Feshbach resonance [12] [13]. While the fast 
sweep regime was analyzed in Refs. [HI [15], the com- 
plementary case of slow sweeping rates, equivalent to a 
many-body LZ problem |16l 117] , has not yet been under- 
stood. 

Having the above setup of two level systems coupled 
to a cavity mode in mind, we consider the Hamiltonian 

\+ N N 
H = At St S+ 1 £ o-* + J* £ (St (j~7 + b at) , (3) 

i=i i=i 

where St creates a photon mode, and of are raising 
and lowering operators of the i— th two level system, 
(cr = (a x ± i(T v )/2, where a x,y,z are Pauli matri- 
ces.) The energy of the photon and the two- level sys- 
tem vary in time as ±At, respectively. The Hamilto- 
nian ([3| is equivalent, up to a gauge transformation, to 

h = -2\t sts + Wo Eilt °t + Eili (& f *r + k + ) ■ 

which represents a generalization of the James-Cumming 
Hamiltonian [T5] to TV two-level systems. Equivalcntly, 
we can think of ^ as an effective Hamiltonian describ- 
ing a Feshbach resonance scenario: representing the spin 
operators in ^ in terms of Anderson pseudospin opera- 
tors, 

of -> - a iX a\ v cr+ -> o^o^, cr" -> &n a iT (4) 

where a\p ah, a^, are the creation and annihilation 
operator for the spin- 1/2 fermions labeled by i, Eq. ^ 
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assumes the form (up to an unimportant constant) 
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(5) 



This is nothing but the Hamiltonian describing the cre- 
ation of molecules out of N fermion pairs in a Feshbach 
resonance experiment |15l 116] (although the single mode 
approximation, i.e. the neglect of bosonic dispersion, 
may be problematic in that case). 

Assuming that the boson level is initially empty, and 
all fcrmions resident in the upper state (on account of 
the energy balance at large negative times), 

(Pb)=0, (<> = 1, i = l,...,N, (6) 

our goal is to compute the asymptotic distribution 

n b (t) = {P(t)b(t)) (7) 

at t — * oo, i.e. the generalization of the LZ transition 
probability P. For N = 1, this task is equivalent to 
the standard LZ problem, whose answer is given by (cf. 

Eq. (|2|) lim^oo nj,(t) = 1 — e However, for TV > 1, 
Eq. pi defines a genuine many-body problem and the 
solution of the corresponding Schrodinger equation be- 
comes progressively more difficult. 

While we do not know how to handle the problem for 
arbitrary N, an approximate solution valid in the limit 
of large particle numbers can be found. At large N, the 
number of produced bosons turns out to be reasonably 
well approximated by 
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lim nb{t) — — 1, 

t— >oo 


< N, 


(8) 


hm n b (t) ~ — ^ — , 


e x ^ iV, 


(9) 






lim n b (t) -> N, 

t — >oo 


e~ > N. 


(10) 



According to these equations, the adiabatic ground 
state (n b t: ^> N) is only realized if A <C irg 2 / log (AT), 
a criterion which is progressively more difficult to satisfy 
as N becomes larger (see Ref. [T3] for a general discus- 
sion of the applicability of the adiabatic limit in large 
systems). This is in marked contrast to the few body 
case, where adiabatic ground state occupancy is granted 

7T0 2 

for large values of the LZ parameter . The observa- 
tion of this difference, obtained for a specific model but 
likely indicative of a more general phenomenon, repre- 
sents the main result of the paper. 

Technically, Eqs. ^ and Eq. ^ obtain by integration 
of a rate equation for the boson occupation number. De- 
noting the occupation of individual fermion states by n/, 



the latter reads 

d t n b = 2ng 2 6(2\t) (n) (1 + n b ) - n b (1 - n f ) 2 ^ , 
n b + Nn f = N 1 (11) 

where the factor S(2\t) accounts the energy balance 
in particle conversion processes, the first (second) term 
on the right hand side describes the creation (destruc- 
tion) of bosonic particles by destruction (creation) of 
two fermions, and the second line enforces particle num- 
ber conservation. Postponing the derivation of Eq. ( |TT| ), 
and the discussion of its range of validity to below, we 
note that upon introduction of a variable 9, such that 
0'(t) — 8(t), Eq. (Ill assumes the form 



dgn b 



■Kg 



(n}(l + n b )-n b (l-n f ) 2 y (12) 



At 9 = 1 (which corresponds to t — > oo) the solution of 
this equation (with boundary condition n b = at 8 = 0) 
reads 



N 



n b 



(.*-.) 



2e"s- 



N 



(13) 



where terms of 0(N 1 ) have been ignored. Taking 
the limit N — > oo at fixed e* g2 f x gives Eq. while 
e^A ^ iv leads to Eq. Although that latter limit 
is beyond the scope of the large N expansion, ^ turns 
out to provide a reasonable (if uncontrolled) approxima- 
tion to n b . 



To actually derive Eq. (Ill, we apply the Keldysh for- 
malism. Defining t = (ti + ^)/2 and r = t\ — t%, we 
denote by G R < A < K {t x , t 2 ) = G R ^ K {t, r) the retarded, ad- 
vanced, and Keldysh fermionic propagators, respectively. 
(For the general definition of Keldysh Green functions 
and notation conventions we refer to the review |20j.) 
Initially all the N fermionic levels are occupied; this cor- 
responds to the bare (noninteracting) propagators 



G«' A (t,r) 



le 



(14) 



where the upper or lower sign in ± and =p are chosen 
for retarded and advanced propagators respectively. We 
aim to compute the boson's Keldysh propagator D K (t, r) 
which, when evaluated at t = +oo, gives the number of 
produced bosons. Initially, however, the boson level was 
unoccupied. Thus 



D»> A (t,r) 



-ie 



i\tr 



(15) 



If the self energy S(£, r) of the bosons is known, the 
Keldysh bosonic propagator can be found by solving 
the Dyson equation, (Dq 1 — S) o D = 1 where Do (D) 
is the bare (dressed) bosonic propagator. Introduc- 
ing the bosonic distribution matrix F(t, t') through [5D] 
D K = D R o F — F o D A , where (A o B)(t x ,t%) = 
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f™ oo dt 3 A(t 1 ,t 3 )B(t 3 ,t 2 ), and D R > A ' K are the retarded, 
advanced and Keldysh components of D, the Dyson equa- 
tion for D K translates to a kinetic equation 

[F ° id t + At] = E A ' — (T, r o F — F o £ a ) . (16) 

To approximately solve this equation, we note that only 





b) 



FIG. 1: a) dominant self energy diagram for the bosonic prop- 
agator in the limit N — > oo with g kept fixed. The straight 
lines are fermions, while the wavy lines are bosons, b) non- 
RPA diagram of lesser order in N. 

interaction vertices ~ g 2 /N accompanied by one summa- 
tion over N fermion states survive the limit N — > oo. In 
practice, this means that only the self energy diagram 
depicted in Fig. [I] a) contributes to the boson self energy. 
Processes such as the one shown on Fig [T] b) are frus- 
trated in that the number of fermion summations does 
not compensate for the number of interaction lines. One 
may verify that the same logics excludes any diagram 
other than the one shown on Fig.[l]a). (For a caveat in 
the argument, see below.) 

The diagram shown in Fig. [I] a) translates to 

£^(t,t') = lg 2 J ^ G R > A {t,e-u)G K {t,co) (17) 
E*(t,0 = 4 E / ^G*(t,e-a/)G*(t, W ). 
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on the right hand side we have switched to a Wigner 
representation, 



(18) 



/OO 
dr G(t, t) e lT 
-oo 

Introducing the spectral function 

A f (t,e) = 2lmG R (t,e) (19) 
we obtain an equation for the Wigner transform of F, 
- 2 i-du 



(d t - Xd e ) F(e) = »- J —A f (e - co) A f (w) [1+ 

/ (e - w) / ( W ) - (/ (w) + / (e — u>)) F (e)] . (20) 

Here f(t, e) is the fcrmionic distribution function, and the 
argument t, identically carried by all Wigner functions, is 



suppressed for brevity. In deriving Eq. ( 20 1 we assumed 



that the Wigner transform of products of operators on 
the right hand side (e.g. (E fl o F)(e, T)) can be replaced 
by the product of the Wigner functions (S(e, T)F(e, T)). 



As discussed a few paragraphs further down, this leading 
adiabatic approximation |20j turns out to be exact in 
our case. Also note that Eq. (20 1 was derived without 



specifying whether the fermionic propagators in Fig.[T]a) 
are bare or dressed. 

Noting that in the distant future fermions and bosons 
become effectively uncorrelated and the energy of the lat- 
ter asymptotes to e = —At, our aim is to calculate the 
bosonic distribution function, nb(t) = rib(t, e = —\t). 
To transform Eq. ( |20| into an equation for rib we use 
the general relations rib = (F — l)/2, rif = (1 — /)/2, 
and note that d t rib = (d t — \d e )F/2. Approximating the 
fermion spectral functions by their bare value, A(e) = 
— 27r<5(e — At), we then readily arrive at Eq. ( [TT] ) , where 
all fermionic distribution functions rif(t) = n/(t,0) are 
evaluated at zero energy. 

Let us examine the status of the various approxima- 
tions used in the derivation of the rate equation: In the 



language of diagrammatic perturbation theory, Eq. ( 20 1 



treats the bosonic and fermionic Keldysh components of 
the self energy operators in a self consistent RPA ap- 
proximation (i.e. a scheme wherein all 'crossing' interac- 
tion lines in the bosonic and fermionic self energy are 
neglected, on account of the limit N — > oo.) Notice 
that a naive interpretation of the large N limit would 
suggest to neglect the fermion self energy altogether: in- 
teraction corrections to the fermion propagator do not 
come with a final state summation and are, therefore, 
superficially of 0(N^ 1 ). However, this argument ne- 
glects that in regimes ^ and ( |To| ) above, the bosonic 
distribution function F ~ rib ~ N introduces additional 
N dependence into the theory. (Physically, the macro- 
scopic occupation of the boson level effectively enhances 
the fermion scattering rate.) This mechanism requires 
us to keep the RPA self energy of the fermionic Keldysh 
Green function. However, the self energy corrections to 
the retarded and advanced fermion propagators, which 
are independent of the bosonic distribution function, are 
indeed large N negligible. This latter simplification jus- 
tifies the above approximation of the fermion spectral 
function by its bare value. (For the sake of completeness 
we mention the existence of non-KPA diagrams in which 
a nominally low power in A^ 1 competes with factors rib- 
[This happens, e.g., in the Keldysh sector of the diagram 
shown in Fig. [l] b) .] These processes are not captured 
in our present analysis which means that the theory be- 
comes effectively uncontrolled once rib ~ N.) 

We finally comment on the status of the leading order 
Moyal expansion (S o F)(e,T) ~ S(e, T)F(e, T) used in 
the derivation. The temporal singularity ~ 5(t) of the 
collision integral makes one worry that this replacement 
may, indeed, not be entirely innocent. While we cannot 
really justify the approximation in the resonant time win- 
dow t ~ 0, we have checked that it does yield the correct 
long time asymptotics ^ when applied to the standard 
LZ evolution equation. 
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It is instructive to reconsider the derivation of Eq. Q 
from a somewhat different perspective: the fact that the 
Hamiltonian (|3| contains the Pauli matrices only in cer- 
tain linear combinations enables us to attack the problem 
by spin algebraic methods. We define an SU(2) algebra 
of spin operators {S z , S + , S~} acting in an N/2 dimen- 
sional Hilbert space as S z = \ Y^i=i a ti ^ ± = J2iLi G± ■ 
Eq. ^ enforces full initial polarization, (S z ) — N/2. 

Since the total number of bosons produced is much less 
than N (the defining criterion of the regime Eq. p|), 
the spin will not deviate much from the vertical direc- 
tion, and it is convenient to employ a Holstcin-Primakoff 
representation: replacing [2T] S + — ► VnIhp, S~ — > 



Nb ( HP , S z 



N/2 



where V HP and b P p 



are the creation and annihilation operators of an aux- 
iliary Holstein-Primakoff boson, the large N limit of the 
Hamiltonian Eq. |3]) reduces to the quadratic form 



H = -\t S f S - Xt P Hp b H p + 9 (& f V HP + 



bb 



HP 



(21) 



The solution of the equations of motion of ( 21 ) then leads 
to Eq. (JsJ _ (In a slightly different context, these equa- 
tions have been solved in [22] , where Eq. Q was derived 
for the first time.) However, the above method does not 
appear to be straightforwardly extensible to the regime 
of large transition rates, Eq. Q. 




FIG. 2: The boson production nt as a function of x = 
exp (Tvg 2 /\) /N. Here N = 100 (empty circles) and N = 500 
(filled circles). The straight line represents Eq. |8| and 
the curve is Eq. j9j. The data is obtained by solving the 
Schrodinger equation for Eq. ^ at g = 1, on the interval 
—40 < t < 40, with the small oscillations in the data being 
the artifact of the finite time interval. 

To check the validity of our results we have run a nu- 
merical test. The above spin representation shows that 
the Hilbert space of the problem is of dimension N + 1 
(much lower than the 0(2^) naively suggested by the 
representation ([3])); this makes a numerical solution of 
the Schrodinger equation feasible. Fig. [2] shows rib as 
a function of x = exp(7r<? 2 /A) /N for N = 100 and 
N = 500. At N — 100 the data is in general agree- 
ment with Eqs. (JsJ) and (pj, at larger N we observe grad- 
ual downward deviations. Preliminary results based on 



a combination of semiclassical ideas and numerical in- 
tegration indeed suggest the existence of corrections in 
ln(A^) (at fixed value x > 1). However, in all our simula- 
tions, the fraction n^/N converged to values below that 
predicted by Eq. (J9|, i.e. our principal observation of 
incomplete ground state occupation remains valid. 

To conclude, we have studied a genuine many body 
generalization of the Landau-Zener problem. Unlike with 
the single particle case, the system does not settle in its 
many particle ground state and a finite fraction of parti- 
cles remains in energetically high-lying states. This phe- 
nomenon makes many body Landau-Zener physics pro- 
foundly different from the few body case. 
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